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基于 核 学 习 方 法 的 短 时 交通 流量 预测 


王 秋 和 莉 ， 李 军 
(兰州 交通 大 学 自动 化 与 电气 工程 学 院 ， 兰 州 730070) 


摘 ww 提出 一 类 基于 核 学习 方 法 的 预测 模型 。 核 递 扒 
二 乘 (KRLS) 基 于 近似 线性 依赖 (approximate linear dependence，ALD) 技术 可 降低 计算 复杂 度 及 存储 量 ， 是 一 种 在 

线 核 学 习 方 法 ， 适 用 于 较 大 规模 数据 集 的 学 习 ; 核 偏 最 小 二 乘 (KPLS) 方 法 将 输入 变量 投影 在 潜在 变量 上 ， 利 用 输入 与 

输出 变量 之 间 的 协 方差 信息 提取 潜在 特征 ; 核 极 限 学 习 机 (KELMJ) 方 法 用 核 函 数 表 示 未 知 的 隐 人 金 层 非 线 性 特征 映射 ， 通 
过 正则 化 最 小 二 乘 算 法 计算 网 络 的 输出 权 值 ， 能 以 极 快 的 学 习 速 度 获得 良好 的 推广 性 。 为 验证 所 提 方 法 的 有 效 性 ， 将 

KELM、KPLS、ALD-KRLS 用 于 不 同 实测 交通 流 数 据 中 ， 在 同等 条 件 下 ， 与 现 有 方法 进行 比较 。 实 验 结果 表明 ， 不 同 

核 学 习 方 法 的 预测 精度 和 训练 速度 均 有 提高 ， 体 现 了 核 学 习 方法 在 短 时 交通 流量 预测 中 的 应 用 潜力 。 
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Short-term traffic flow forecasting based on kernel learning methods 


Wang Qiuli, Li Jun 
(School of Automation & Electrical Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China) 


Abstract: Based on the powerful nonlinear mapping ability of kernel learning, proposing a class of kernel learning method for 
the short-term traffic flow forecasting. Kernel recursive least squares (KRLS) method using Approximate Linear Dependence 
(ALD) technique can reduce the computational complexity and storage capacity, the KRLS method is an online kernel learning 
method and is suitable for training on large-scale data sets. Kernel partial least square (KPLS) method utilizes the covariance 
between input and output variables to extract latent features. Kernel extreme learning machine(KELM) method uses the kernel 
function to substitute for the unknown nonlinear feature mapping of the hidden layer, in addition, the output weights of the 
networks can also be analytically determined by using regularization least square algorithm, hence KELM method provides 
better generalization performance at a much faster learning speed. In order to verify the validity of the proposed kernel 
learning method, the employed ALD-KRLS, KPLS and KELM methods were respectively applied to different traffic flow 
forecasting instances in different area, compared to the other methods under the same conditions. Experimental results show 
that the proposed kernel-based learning methods have higher forecasting accuracy and have improved training speed in the 
short-term traffic flow forecasting. 


Key Words: kernel learning methods; short-term traffic flow:; forecasting 


0 引 器 学 习 方 法 近年 来 在 短 时 交通 流 预 测 方面 已 取得 许多 成 功 应 
用 。 文 献 [3~5] 将 神经 网 络 应 用 于 交通 流 预测 模型 。 文 献 [6~8] 

交通 流量 预测 作为 智能 交通 系统 (intelligent transportation 分别 将 SVM、 最 小 二 乘 支持 向 量 机 (LSSVM) 应 用 于 短 时 交通 流 
systems, ITS ) 中 的 重要 组 成 部 分 ， 在 城市 路 网 交通 信号 控制 、 单 步 预测 中 。 极 限 学 习 机 (ELM) 外 随机 选择 SLFNs 的 隐 含 层 节 
拥塞 控制 、 驾 驶 路 线 规划 等 方面 发 挥 重要 作用 03。 交 通 流量 预 。 点 及 相应 参数 ， 具 有 极 快 的 学 习 速 度 ， 当 ELM 的 隐 含 层 节点 
测 可 为 旅行 者 提供 实时 交通 信息 以 节约 出 行 时 间 与 成 本 ， 更 有 未 知 时 ， 采 用 核 矩 阵 方法 替代 之 ， 形 成 KELMIUI 方 法 ， 文 献 
助 于 交通 决策 部 门 制定 积极 主动 的 交通 管理 策略 ， 以 提高 城市 。 [10~13] 成 功 将 ELM、KELM 应 用 于 时 间 序 列 预测 中 。 
交通 路 网 的 运行 效率 。 考虑 特征 提取 技术 ， 文 献 [14，15] 分 别 将 SVM 与 主 元 分 
神经 网 络 与 支持 向 量 机 (support vector machines，SVM) 等 析 、 核 主 元 分 析 结 合 的 PCA-SVM、KPCA-SVM 方法 应 用 于 短 
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时 交通 流 预测 中 ， 在 一 定 程度 上 均 有 效 提高 了 预测 精度 。 
文献 [16] 给 出 了 KPLS 方法 ， 它 是 线性 PLS 方法 在 非 线 性 
特征 空间 上 的 延伸 。 考 虑 在 线形 式 的 核 学 习 方法 ， 文 献 [17] 给 
出 了 KRLS 方法 ， 它 基于 稀疏 技术 限制 核 矩 阵 的 规模 ， 以 适应 
较 大 数据 集 的 训练 。 文 献 [18] 将 KRLS 方法 用 于 时 间 序 列 预 测 
中 ， 具 有 较 好 的 预测 效果 。 
鉴于 SVM、ELM 等 方法 在 短 时 交通 流 预测 等 方面 的 成 功 
应 用 , 考虑 到 KELM 等 核 学 习 方 法 在 短 时 交通 流 预 测 领 域 的 研 
究 成 果 相 对 少见 ,针对 上 述 问 题 ,本文 提 出 基于 KELM、KPLS、 
KRLS 的 一 类 核 学 习 预 测 方法 ， 以 进一步 提高 预测 精度 及 提高 
计算 效率 。 
为 验证 所 提 基 于 KRLS、KELM、KPLS 方法 的 预测 模型 的 
有 效 性 ， 将 它们 应 用 于 某 地 区 的 实测 数据 实例 中 ， 在 同等 条 件 
下 ， 不 仅 对 所 给 出 的 三 种 不 同 核 学 习 方 法 进行 比较 ， 而 且 还 将 


要 次 莉 ， 等 : 学 习 方法 的 短 时 交通 流量 预 六 


Pr yr 
p=| : Y= (5) 
Bs J; yw ys 
其 中 : 为 SLFNs 的 隐 含 层 输出 矩阵 ， 玖 的 第 i 列表 示 与 输 
入 xz 相关 的 第 i 个 隐 含 层 节 点 的 输出 向 量 ; 第 i 行 表 


示 与 输入 x, 相关 的 隐 含 层 特征 映射 。 


若 激活 函数 G(x,;w,,b) 在 任意 区 间 上 无 限 可 微 , 则 当 隐 含 
层 节 点 数目 地 小 于 样本 N 时 ，SLFNs 仍 能 以 极 小 的 训练 误差 
逼近 训练 样本 外。 此 时 托 阵 如 并 非 方 阵 , 因此 存在 记 ， 放 ,和 忆 
使 得 
疏 Ge Wb. Bb) —Y 
(0) 
= min HOw wi bs be)p 7)| 


忆 此 ， 与 一 般 的 前 馈 神 经 网 络 方法 调整 各 层 网 络 权 值 的 训 


与 现 有 SVM 等 预测 方法 进行 实验 对 比 ， 以 进一步 拓宽 核 学 习 
方法 在 智能 交通 系统 中 的 应 用 潜力 。 
1 ” 核 学 习 方法 

给 定 训练 数据 {x,,y,}e R”"xR"，i=1,2,.…,N ， 构 建 秆 阵 
和 ER 了 <Rww 。 将 数据 x 映射 至 高 维特 征 空间 中 ， 即 
:xeR" ->G(x)eFcR”， 为 避免 特征 空间 上 高 维 计算 的 困 
难 ,考虑 定义 满足 Mercer 条 件 的 核 函 数 K(x,,x,) = G(x,) G(x))， 
则 有 以 k(x,,x,) 为 元 素 构成 的 核 矩 阵 天 -GrEG ， 其 中 
B=[Gx) GX,)...Gxn)] 。 
1.1 KELM 方法 

KELM 是 ELM 方法 在 非 线 性 特征 空间 上 的 延伸 , 首先 考 
虑 对 于 含 启 个 隐 舍 层 节 点 的 SLFNs， 其 网 络 输出 为 
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y(x)= yp = ypGCcwpDix ER"”",B,sR"(D) 
过 1i=1 


其 中 : p 表示 连接 第 i 个 隐 含 层 节点 与 输出 层 节点 的 权 值 向 量 。 
凡 G9 表示 第 i 个 隐 含 层 节点 的 输出 函数 G(x;w,,b) ， 由 隐 含 层 
节点 的 类 型 确定 。 若 隐 含 层 节点 为 RBF 节点 ,其 激活 函数 形式 
汉 


G(xX;w,,b,) =exp(-b -wl ) Q2) 


若 SLFNs 能 以 零 误 差 通 近 N 个 训练 样本 ， 即 


jo -|=0， 则 存在 p,，w, 和 5b 有 如 下 表达 式 : 


BG :wb) = = (3) 
矩阵 形式 可 表示 为 
HP=Y (4) 
h(x,) G(wi:xi+b) G(wi :xi+b;) 
H=| : |= : 四 : 
h(x,) G(wi-xw+b) 的 让 G(wi xXw+b;) NxN 


练 学 习 不 同 , ELM 方法 仅 需 训练 网 络 的 输出 权 值 矩阵 8 ， 其 学 
习 过 程 相当 于 求解 式 (4) 的 方程 组 Hp = 了 的 最 小 二 乘 解 名 ， 即 
有 


p=HY 0) 
为 进一步 提高 网 络 的 泛 化 能 力 及 和 矩阵 求解 的 稳定 性 ， 根 据 
岭 回 归 思 想 引 入 正则 化 参数 ， 对 p 的 解 可 进行 一 定 约束 。 医 


此 ， 式 (7) 的 求解 可 进一步 表示 为 : 
IY 
-pT ,i (8) 
p=H [an + 了 


相应 地 ，ELM 网 络 的 输出 为 

y(x)=h(x)p (9) 

在 式 (8)(9) 的 求解 过 程 中 , 进一步 考虑 用 核 函 数 表示 未 知 的 

隐 含 层 非 线 性 特征 映射 ， 引入 满足 Mercer 条 件 的 核 函 数 ,定义 

矩阵 天 = BET, 其 元 素 为 [区 = K(x,x))=h(%)ohX)) 。 此 时 ， 
KELM 网 络 输出 为 


区 


-1 
y(x)=h(x)p -=n (Ha + 5] 了 


风骨 所 (10) 
= (x Y 
a 


k(xX, XN) 


1.2 KPLS 方法 

KPLS 将 高 维特 征 空间 看 做 原始 空间 的 对 偶 ， 基 于 核 变 换 
技术 在 对 偶 空间 中 建立 线性 PLS 回归 模型 。PLS 方法 利用 输 
入 和 和 输出 y 矩阵 的 协 方差 选择 潜在 特征 变量 , 其 回归 模型 可 
写成 矩阵 形式 : 


Y=XB+F (11) 
其 中 : B 是 维 数 为 mxn 的 回归 系数 矩 阵 ，p 为 维 数 为 N xn 的 
残 差 矩 阵 。 
文献 [19] 将 线性 PLS 方法 延伸 至 高 维 非 线 性 特征 空间 ， 
成 KPLS 方法 ， 其 算法 实现 如 下 : 
a) 随 机 初始 化 输出 矩阵 y 的 湾 在 向 量 u ， 并 定义 核 矩 阵 
玉 =GIG 。 


NY 


录用 入 
b) 计 算 输 入 矩阵 这 的 潜在 向 量 4 ,并 进行 归 一 化 处 理 , 即 ; 
1 GTOx = Ku, 全 (12) 
co) 计 算 载荷 向 量 c 。 
c = 了 It (13) 
gd) 计算 x 并 进行 归 一 化 处 理 ， 即 : 
u=Yc, u 二 (14) 
lal 


6) 重 复 b)~e)， 直 至 收敛 。 
用 对 核 矩 阵 K 及 输出 矩阵 


进行 收缩 映射 ， 即 


K=(I-tt" )K(I tt )=K ttIK — Ktt™ +tt'Ktt’ (15) 


Y=Y-tt'Y (10) 
可 归 系 数 矩 阵 如 = BU(TTKU)'TTY 。 因 此 ， 多 维 输入 、 
多 维 输出 的 KPLS 回归 模型 在 训练 数据 集 上 的 预测 输出 为 
Y=GB= KUCTKU TY =TT'Y (17) 
其 中 :矩阵 TU 的 各 列 由 向 量 w 构成， 矩阵 TT 的 各 列 由 向 量 1 构 


为 消除 回归 模型 中 偏 置 项 的 影响 ,还 应 对 映射 到 特征 空间 
中 的 数据 进行 中 心 化 处 理 ， 需 要 对 训练 数据 及 测试 数据 所 构成 
的 核 矩 阵 分 别 进行 如 下 变换 ， 即 : 


1 1 

R=- (1 g(r | (18) 
> 让 1 .二 (19 
K,= K, -Tylak, 工交 ll ) 


其 中 : K, 为 未 进行 中 心 化 的 测试 数据 核 矩 阵 , 维 数 为 N, x N,N, 
为 测试 数据 集 的 数据 大 小 , 1 是 元 素 为 1, 长 度 为 N 的 列 向 量 
J1x 是 元 素 为 1， 长 度 为 N 的 列 向 量 ，7 为 单位 阵 。 
考虑 用 于 时 间 序 列 预测 建 模 的 单 输出 情形 ， 此 时 回归 系数 
向 量 4=UCTKU)-TI7 ， 其 元 素 为 di=1..,N 。 则 KPLS 预 
测 模型 的 式 (17) 可 进一步 表示 为 


oe Ya (20) 


式 (17)， 若 仅 提取 前 p 个 潜在 变量 ， 在 非 线性 特征 空间 
中 ， 将 KPLS 预测 模型 解释 为 线性 PLS 回归 模型 ， 即 
y(z)=cn(z)+…+cop(z)=ct(z) (21) 


pp 


其 中 : {i (x))*, 是 x 的 前 p 个 潜在 变量 ， 即 非 线性 主 元 上 的 投 


1.3 KRLS 方法 
在 训练 阶段 ， 假 设 获取 1 时 刻 之 前 的 训练 数据 
由 训练 数据 的 子 集 可 构成 字 


人 后 和» 
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Ean! 


-{E j5 ， 其 中 元 是 所 选择 的 训练 数据 ， s， 是 字典 集 全 


Is 


的 数目 。 特 征 空间 中 ，{g(%,)} 忆 是 线性 独立 的 数据 向 量 。 当 新 
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数据 x 出 现时 ， 判 断 g(x,) 是 否 满足 ALD 条 件 ， 若 5 >v， 则 
将 x 加 入 字 与 万 中 。 

为 考虑 是 否 将 数据 x, 加 入 字典 集合 中 ， 需 求 出 满足 ALD 
条 件 的 系数 向 量 a =(a,,.…,a, ,)"， 即 有 


Sap) -G(s) 


其 中 净值 参数 k(x,,x),) = G(x,)™G(x,) 的 取 值 大 小 影响 稀疏 化 程 
度 。 考 处 k(x,,x)) = x,) G(x,) ， 由 式 (22) 得 到 


5 -nin Saka)- 


i,j=1 


底 


0 = min <v (22) 


1 


(23) 


2 a (zx)+ et 
癌 


=min {a'R, a —2ak,, (x, ) 十 后 


i, 7 =1,...,S, 1 
求解 式 (23)， 可 得 到 最 优 系数 向 量 a 及 ALD 条 件 : 
a, = Kk _ (x,) ， 
6 =k, k(x) a <v (24) 
党 助 ALD 条 件 ， 直 至 上 的 所 有 数据 均 能 由 字典 集合 中 
数据 向 量 的 近似 线性 组 合 来 表示 。 由 式 (22) 可 推 得 
gr) 人 rs 9 


其 中 如" 为 残 差 成 分 向 量 ，a,, 是 系数 向 量 a 的 第 7 个 元 素 。 延 
伸 至 特征 空间 ， 定 义 矩 阵 @=[G(%) GE,)...G(£,)] 
Ge =[G(x*) G(x*)...G(x*)] 。 式 (25) 的 矩阵 表达 为 


B=-BAT+HOD™ (26) 
当 y 足够 小 时 , 式 (26) 的 近似 表达 式 为 = A,KKAT， 其 中 


[K,];; = K(xi, xX,)), [A] = 4» b,j=1,...,to 


“， 由 RLS 算法 


Dw—Yy 


1 1 


在 时 刻 L(w) = 了 (f(x)-)= 


可 知 ， 最 小 化 误差 平方 和 满足 下 式 : 


zw)=y(7(z)-y)=|em-7 


1 


pl 
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其 中 矩 阵子 =(y4;.…;y,) ， 维 数 为 1xn。w 为 权 值 系数 向 量 。 | 


表示 2 范 数 。 由 式 (27) 得 w 的 最 优 解 为 


w, = Yag(x,) =Gu (28) 
其 中 w=(w，.,a)7 ， 由 式 (26)(C27) 得 到 
raw=lka-z 9 


式 (29) 进 行 稀疏 化 处 理 ， 令 字典 数据 向 量 形成 的 核算 阵 开 
逼近 天 ， 开 是 由 字典 集合 中 的 所 有 8 个 元 素 构成 的 核 矩 阵 ， 


借助 ALD 条 件 从 训练 数据 中 构建 字典 集合 。 
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式 (26) 代 入 式 (28)， 则 w=@Ba =GArwu =GBa ， 其 中 四 
j= Ara 为 st 个 “ 约 简 ”系数 构成 的 向 量 。 因 此 ， 式 (29) 变 形 


it 
时 RK- (P,, -pa 0 PP， )(ALY,, +ag,) 


类 机 本 38 
=C + 开 (Pad -BarRC 9 
ro) -leoa-rl rT 60 -6 tk pg -E(w) 6) 

式 (30) 的 最 小 化 求解 为 输出 字典 ，& 。 
(AR)Y =R- (ATA) AT 羡 gj 迁 代 计算 Step 2~ Step 4， 直 至 所 有 训练 数据 依次 完成 。 
0 和 测试 阶段 ， 对 已 训练 的 KRLS 模型 进行 测试 ， 其 预测 输 


多 维 输入 -输出 ALD-KRLS 方法 的 具体 实现 如 下 : 出 如 式 (39) 所 示 。 
引 训 练 阶段 ,给 定 阔 值 参数 ， .时 刻 寺 1 时 ,获取 数据 (x,,y)， 


3 =k, (x, ) Ci G39) 
K=[k,]; 天 =[1/k,]; a =(% /hk) P=[1]; $=]。 


若 训 练 数 据 的 数目 为 N， 则 KRLS 算法 的 计算 复杂 度 为 
b) 令 六 寺 1， 得 到 新 数据 (x,,y,) ， 计 算 OCV2) 。 


-TT 


En) | 局 击 志 丰 天 雹 二 5 机 2 G32 2 短期 交通 流 预 测 实例 


0O) 计 算式 (24), 如 果 5 >v ,将 x, 加 入 字典 中 ,字典 矩阵 EK, 按照 时 间 序 列 建 模 的 方式 ， 将 本 文 的 不 同 核 学 习 方法 应 用 
也 相应 增长 。 即 有 忆 =D_U{E}， s=s+l1; 更 新 计算 展 - ; 于 短期 交通 流 预测 实例 中 ， 所 建立 的 预测 模型 为 
K, k(x,) P(t+h)=f (x,) (40) 
k(x, ) k, 全 G3) 其 中 :h 表示 预测 步 长 , x 是 历史 时 间 序 列 值 构 成 的 多 维 输入 向 
EK, _ 工 | 5 天 1 TaaT a, 量 x =[y%.…y4] 为 嵌入 维 数 。 对 于 单 变量 交通 流 时 间 序列 ， 
0 -a 1 鸭 建 训练 数据 对 {x,y} eR"xR, i=1,2..,N, 此 时 模型 输出 
于 x, 是 字典 中 的 数据 向 量 ， 相 应 有 : 维 数 n=1，/(-) 则 分 别 对 应 于 KRLS、KPLS、KELM 方法 ,并 
了 在 同等 条 件 下 ， 与 KPCA-SVM、SVM 等 方法 对 比 。 
| Da 1 (34) 人 与 方法 对 
0 1 1 实验 采用 正则 化 均 方 根 误差 aormalized mean square root 
P 更 新 为 error, NRMSE) 以 及 均 方 根 误差 (mean square root error, RMSE)、 
Pp- (ATA _ 区 "| G35) 百分比 误差 (mean absolute percentage error, MAPE) 作 
0 1 为 衡量 指标 。NRMSE 表达 式 如 下 : 
其 中 0 是 零 向 量 。 % 
式 (24) 得 到 or 天 ,= 无 ,(z)7， 因 此 人 更 新 为 >》 (0 一 幼 ? 
CE NRMSE_」 与 (41) 
6 =R'(A'A,) AZ (N, —D var(y) 
gil (MA) Ds w 
=K, 其 中 : $ 表示 模型 的 预测 输出 ，y 为 实际 输出 ， var(y) 表示 序列 
(36) 方差 ， N 为 测试 数据 数目 。 实 验 中 的 高 斯 核 函数 形式 如 下 : 
CQ _ gq, 一 大- (x, ) Ci 1 
了 | 5 | k(xix))=exp{-|s; -x|/20°} (42) 
1 ~ i 
二 (4 —k,_ (x,) 5 i 2 
0 其 中 : c 为 高 斯 核 函数 的 宽度 。 
输出 字典 Dp ，&。 2.1 实例 1 
四 如 采 5<py，D=D,，s=s， 久 =K，,。 更 新 选取 北京 某 公 路 观测 站 的 实测 交通 流 数 据 ， 观 测 总 时 段 为 
T A SN 4 天 共 96 小 时 ， 以 15 4 为 采样 间隔 ， 记 和 384 组 数据 。 
二 二 ] 。 按 箱 阵 求 逆 引 理 更 新 : 天 小 时 ， 以 15 分 钟 为 采样 间隔 ， 记 录 了 组 数据 
实验 中 取 前 285 组 数据 为 训练 数据 ,后 92 组 数据 为 测试 数据 ， 
元 _PaaP, G7 选取 先入 维 数 m=4。 由 交叉 验证 法 确定 不 同 核 学 习 方 法 的 预测 
ltaPna, 模型 参数 。KELM ，KPLS， KRLS，KPCA 均 各 自选 取 高 斯 
核 函数 ， 核 函数 宽度 oc-1 。KELM 正则 化 参数 y=50 ; KPLS 
着 定义 局 让 全 全， 则 局 = Ra ， 下 Cn)= 关 w 更 中 汪 丰 全 县 下 日， 
1+ pa 中 潜在 变量 数目 p=15 ; ALD-KRLS 最 大 字典 容量 smax=200， 


新 计算 & 如 下 : 阔 值 yv=0.1 。 对 比 实 验 中 ，SVM 方法 使 用 Libsvm 软件 完成 ， 
惩罚 因子 C=0.5，g =0.01; KPCA-SVM 方法 中 ， 选 取 非 线 


录用 稿 
性 主 元 数目 为 15，SVM 的 核 函 数 为 线性 核 函 数 。 结 果 评 价 1 
表 1、2 列 出 。 
表 1 不 同 核 学 习 方法 的 提前 15 min 的 预测 结果 对 比 
预测 方法 MAPE RMSE NRMSE 
SVM 33.8086 23.4548 0.2693 
KPCA-SVM 33.9102 23.3346 0.2679 
KELM 33.4629 23.2607 0.2677 
KPLS 33.4437 22.5822 0.2593 
KRLS 33.8606 22.5210 0.2506 


表 2 不 同 核 学 习 方 法 的 提前 30 min 的 预测 结果 对 比 


预测 方法 MAPE RMSE NRMSE 
SVM 43.4466 30.4298 0.3476 
KPCA-SVM 43.9819 30.0215 0.343 0 
KELM 45.2888 29.8679 0.3412 
KPLS 47.1310 29.2825 0.3345 
KRLS 42.0246 28.2343 0.314 4 


测 中 均 取得 了 较 好 的 预测 效果 ， 
图 1、2 分 别 给 出 了 不 同 核 学 习 方 法 在 测试 集 上 进行 提前 


15 min、30 min 预 涡 


E 确 地 预测 实际 交通 流量 


KRLS 预测 效果 最 佳 。 


的 对 比 曲 线 。 


多 


表 1、2 可 看 出 ， KELM、KPLS、KRLS 在 交通 流量 


时 了 预 


1、2 看 出 ， 不同 核 学 习 


呈现 出 好 的 预测 效果 。 图 
在 训练 数据 集 上 其 均 方 误差 MSE 的 变化 曲线 ， 


方法 具有 较 快 的 收敛 速 度 。 
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图 1 基于 核 学 习 方 法 的 提前 15 min 的 预测 结果 
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图 2 基于 核 学 习 方法 的 提前 30 min 的 预测 结果 


值 ,预测 误差 波动 相对 较 小 ， 
3 还 给 出 KRLS 方法 进行 单 步 预测 时 ， 
由 图 可 知 , KRLS 
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对 
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2.2 实例 2 


150 200 


样本 


训练 集 上 KRLS 方法 的 收敛 曲线 


50 100 


实例 选取 英国 某 地 区 交通 局 提供 的 交通 流量 数据 ， 采 样 间 
隔 为 15min,， 采样 时 间 为 2011 年 3 月 ， 本 实验 截 选 了 336 个 时 


间 点 。 选 取 骨 入 维 


数 m=3。 实 验 中 取 前 235 个 数据 作为 训练 集 


输入 ， 最 后 100 个 数据 作为 测试 样本 。 


参数 设 定 如 下 所 示 。KELM ，KPLS， 
SVM 均 各 自选 取 高 斯 
化 参数 y=100; KPLS 中 潜在 变量 数目 p=20; ALD-KRLS 
大 字典 容量 smax=90， 阐 值 y=0.1。 对 比 实验 中 , SVM 方法 惩罚 
因子 C=0.5，g =0.01; KPCA-SVM 方法 中 ，KPCA 的 非 线性 


KRLS, KPCA, 


主 元 数目 为 10，SVM 选取 线性 核 函数 。 

表 3 不 同 核 学 习 方法 的 提前 15 min 的 预测 结果 对 比 

预测 方法 MAPE RMSE NRMSE 

SVM 16.4543 66.5587 0.698 7 

KPCA-SVM 15.9120 62.6299 0.6261 

KELM 15.698 1 61.4080 0.6452 

KPLS 15.3897 58.4583 0.5878 

KRLS 13.8856 55.0122 0.4920 


表 4 不 同 核 学 习 方法 的 提前 30 min 的 预测 结果 对 比 


预测 方法 MAPE RMSE NRMSE 
SVM 20.5209 76.0844 0.7704 
KPCA-SVM 19.665 1 74.9778 0.6839 
KELM 18.9895 72.1687 0.6648 
KPLS 17.1630 69.6158 0.6210 
KRLS 15.652 1 65.4789 0.5482 


图 4.5 分 别 为 不 同 核 学 习 方 法 分 别 进行 提前 15 min、30 min 


的 单 步 、 多 步 预 测 时 在 测试 集 上 的 预测 结果 对 比 曲 线 , 由 图 
本 文 的 核 学 习 方 法 呈现 出 较 好 的 预测 效果 。 
, 不 同 核 学 习 方 法 分 别 进行 提前 15 min、30 min 


5 的 结果 看 出 ， 
在 测试 身 


uy 


4、 


预测 的 预测 结果 评价 指标 对 比 则 由 表 3、 表 4 给 出 。 从 表 3、 表 
4 可 以 看 出 , 与 SVM 方法 相 比 , KELM 方法 的 预测 精度 略 有 提 
升 ，KPLS 方法 次 之 ，KRLS 方法 则 显示 出 了 较 高 的 预测 精度 。 


通过 算法 训练 的 收敛 时 间 对 比 可 知 ,，KELM 与 KRLS 方法 的 算 


法 收敛 速度 较 快 ; KPLS 预测 方法 的 训练 速度 较 慢 ; ALD-KRLS 
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方法 具有 自 适 应 在 线 训练 特性 , 其 收敛 速度 快 ,预测 精度 最 高 。 
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图 4 ”基于 核 学 习 方 法 的 提前 15 min 的 预测 结果 
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图 5 基于 核 学 习 方 法 的 提前 30 min 的 预测 结 
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针对 交通 流 短 期 预测 , 分 别提 出 基于 KELM、 KPLS、KRLS 
的 核 学 习 预 测 模 型 ，KELM 方法 以 核 映 射 代替 了 传统 ELM 的 
随机 映射 ， 提 高 了 网 络 的 泛 化 能 力 和 学 习 速 度 。KPLS 方法 
有 特征 提取 能 力 强 、 算 法 实现 容易 等 优点 。KRLS 方法 具有 时 
变 自 适 应 学 习 特性 ， 通 过 稀 跑 化 算法 ， 可 控制 核 箱 阵 的 大 小 ， 
从 而 提高 计算 效率 。 实 验 结果 表明 ， 将 不 同 核 学 习 方法 应 用 于 
短期 交通 流 预测 中 是 一 种 可 行 且 有 效 地 选择 , 另 一 方面 ,KELM 
与 KPLS 方法 均 适宜 于 离线 训练 ， 其 不 足 在 于 当 训 练 数 据 集 的 
数据 超过 一 定 规模 时 ， 会 引起 核 矩 阵 的 运算 困难 。KRLS 方法 
则 适宜 于 在 线 训练 ， 使 其 能 用 于 较 大 规模 数据 集 的 训练 。 
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